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Abstract 

We prove that it is possible to approximate Z by an integral G with 
respect to the standard Gaussian measure such that the approximation 
error \G — Z\ = 0(N~ 1 ) is small, and the value of the integral 



with D = (1 — iB/y/N), is finite and oscillates rapidly in A and 
B. Straightforward summation of Z is impossible because over the 
huge interval n G [N — y / 2iVTog _ /V, N + \/2N log iVJ , the summands 
contribute significantly. We discuss approximation for some oscillating 
series related to Z. 



1 Introduction 

When a mathematician says "consider the asymptotics of a state ip%(x) as 
the Plank constant U goes to zero" , this makes no rigorous sense in physics, 
because h is a constant. The difference between the mathematical infinity 
and huge physical quantities such as the Avogadro number Na or h' 1 is 
that logA a w 23, and log/T 1 w 37 are rather moderate numbers, if they 
are used as multipliers. Therefore, for mathematicians, 1/logA 7 ^ = o(l) 
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and 1 /log ft 1 = o(l) (if they keep Na — > 00 and 1/ft — > 00 in mind), 
and in physics, 1/logA^ = 0(1), 1/logft" 1 = 0(1), log N A /N A = O(N^) 
and ft/log tT 1 = 0(h), e. t. c. Fortunately, a rather general property of 
the nature is that the logarithms of overwhelming physical quantities are 
moderate numbers. On the other hand, the logarithms of large parameters 
occur quite rare in asymptotic expansions. But we will see that crucial 
estimates of the asymptotic expansion of Z relay on the fact that log AT is 
finite. 

The sum Z represents some quantities related to explicit solution of 
a problem on interoferometrical detection of gravitational waves (see pi- 
ll]), where N is the average number of photons generated by a laser, and 
A = A(t), B = B(t) are small time-dependent functions describing interfer- 
ation of photons splitted into two arms of the measuring device. The main 
computational difficulty is that 2y/2N\ogN ~ 10 12 summands centered at 
n = N = 10 23 give a relevant contribution to Z. Thus the straightforward 
computation of Z is impossible. 

In this paper we prove that the series Z(A, B, N) can be approximated 
by the Gaussian integral, because in terms of the new variable x n = ^^jf, 
dx n = ^7= (n — 0, 1, . . .), the Stirling expansion of Gamma-function implies 
the Gaussian approximation of the Poisson probability distribution Pjv(n): 

^^ N i--Aj£^r o(N) ^ (Li) 

if x n /VN is small (see Appendix 3, (A. 6), (A. 7)). 

On the other hand, by the trapezoidal approximation of the integral, it 
is possible to estimate the sum 

Z = £P„(n)e""<»\ S N (n) = ^L+ ^ 



VN 2 v / iV3 

over the interval B N d = [N - ^2N log N, N + ^2N \ogN\ cN. In this case, 
max n£Bjv |a; n | = y/2 log N. More precise, to estimate the sum Z, we must 
justify several approximations: 

- estimate the error of the restriction to Bn of the sum Z; 
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Figure 1: The picture shows the Poisson distribution P/v(n) for N = 10 4 , 
which practically coincides with the Gaussian distribution iV(10 4 , 10 4 ), and 
the approximation error (double-well curve scaled by 10 5 ), representing the 
difference between the Poisson and normal distributions with the Stirling 
corrections. 

- estimate on Bjv the error of approximation of the Poisson distribution 
Pn(P) by the Gaussian one as in Eq. (1.1); 

- estimate the approximation error of the replacing of the sum over 
by the Gaussian integral over the set 

D N = f [-y/2logN, v/21ogiV] c R 

- expand (1.1) in 1/y/N on and estimate the approximation error; 

- extend the Gaussian integral to R and evaluate it analytically. 

First, in Section 2, we cut out the left and right tails of the probability 
distribution P/v(n) outside of Bn, and then in Section 3, we approximate the 
sum by an integral over the interval DjvCR corresponding to B^\ 

Z =7T^ £ (fN(x n )e iS ^ + f N (x n+ i)e iS ^ n+1) ) + Vn = 
2ViV Sn6Bjv V / 

= / dxf N (x)e^ x+ ^wK V N + e N , (1.2) 
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where and r]^ are the approximation errors, <Jn{x) = Sn(xvN + N) and 
/n(x) = \/N Pn(x\/N + N). More precise, for the total approximation error 
VN + £n = Z - J dx f N (x)e iaN ^ x \ we have 

\Vn\ + \cn\ < p N(n) + 

n£]Sl\B N 

+^E{^l4 2) (y)l + 2i^|Pif ) (y)5S ) (y)|+i^|PN(y)5g ) (y)|} > (1.3) 

o ~ I vein yei n yei n ) 

def 

where I n = [n,n + 1] (see Appendix (A.l) and (A. 5)). 

In what follows we prove that |ejv| + \wn\ = 0(N~ 1 ), and Eq. (1.1) de- 
scribe convenient analytical approximation of /n(x) on by the perturbed 
Gaussian distribution (see (A. 6) and (A. 7)). 

This implies that the sum Z can be approximated by the Gaussian integral 

x 2 £6VN 12iV 



r 2 pfi\/Af VAN , 

\ dxe~ x -^ iS N (VN X+ N) + Q{N -i : 



Eventually, the last integral can be extended to R and evaluated explicitly 
up to terms of order 0(iV _1 ): 

i(A+B/2)VN (A+B) 2 / ( A -\- n) 3 \ 

7 - e -2 ( i-iB/VAD ( i _ i SfhZ^t ) + OiN- 1 ). (1.4) 

1-iB/y/N V 6 V^ / 

The last estimate of the error is proved by a generalization of the Komatsu 
inequality (3.3) proved in Appendix 4. This asymptotic expansion (1.4) is 
the main result of the note. It can be applied in many different ways. 

In Section 4, we present similar estimates for several related oscillating 
series. First, it justifies the asymptotic expansion of the sum 

Zp = e — -F ( —= J e 2jv3/ 2 
^ n! vv^v; 

where is the Fourier transform of any finitely supported function F(p) 
such that (diamsuppF+ \B\) 2 < logiV: 

g 2(l-iq/v^V) 1 _ ^ 



1 /■ ~ e HP+ 

V2vr ^fl A _ 



iq/VN ^ " 6^ 



It can also be applied for the correlation and Fourier analysis of Z(t) in 
the case, when A — A(t) and B = B{t). 

Finally, we note that the series Z(x,t,h~ 2 ) d = iph(x,t) represents an 
explicit solution of the Cauchy problem for the Schrodinger equation with a 
special periodic initial condition: 

mi, = = exp{-/T 2 (l - e ix ' n )}. 



2 Estimate of the weight of tails 

The total weight T N of the left and right tails of the probability distribution 
-Pjv(n) on the complement of the interval B N can be estimated by use of the 
expansion (1.1) which is clearly valid on this set. By definition, we have 



T N d ^l- £ P N (n) = l- £ ^=f N {x n ). (2.1) 

neB N x n eD N V jV 



The trapezoidal approximation (see Appendix 1, (A.l) for the estimates of 
the approximation errors) implies the following integral representation of the 
sum: 

X —nc T fN{x n ) = dxf N (x)+~f N , 
. v N Jd n 



where |7jv| < | J2 n eB N max xe/„ \Pn\ x )\. Note that the function Pn(x) has 
just one extremum on the set B^. Hence its first and second derivatives 
may have no more then two and three critical points respectively. Hence the 
function \P^ (x)\ is monotone on all intervals I n with the exception of the 
three (see Fig. 2). 

Therefore, we can apply the estimate (A. 5): 



\ lN \< £ max|p( 2 )(;r)|<3 £ \P®(n)\ = O(N^) 

n£B N Xiz n neB N 

Now the expansion (1.1) ensures the equality: 
/ dx /n(x) = — / dxi 

JD N V27T Jd n 



x 3 x 4 

■ e ~ 



1 J % L _J_ I 

T 2VN T 12V ) 
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Figure 2: The derivatives of the Poisson distribution become more flatten as 
N increases. Each derivation add a new critical point. The picture shows 
Pjp (n) with two critical points, and 100 * Pffi (n) with three critical points, 
for N = 10 4 . 



, _■£_( 3x-x 3 
axe 2 1 =- 



r v +0(N- 1 ) = 1 + 0(N- 1 ) 

because the functions x and x 3 are odd. The estimate of the approximation 
error follows from the Komatsu inequality || (see Appendix 4): 

4e~^(Va 2 + 4 + a] <[ dxe~^ < 4e"^ ( Va 2 + 2 + 

\x\>a V 



For a = a/2 log N, we obtain 
1 



2tt Jr\d 



, ( 3x - x 3 
axe 2 1 



6VN 



1 



4 ^ 
< , e 2 



2tt 



Va 2 + 2 + 



-i 



a=y/2\ogN 



2-K JR\D N 

0{N~ X ) 



dxe 2 < 



These estimates explain the specific choice of the sets D N and B N . Therefore, 
the total weight Tn of the tails on the complement of equals Oi^N^ 1 ): 



T, 



N 



1 - £ P N{n) = 1 - E -nxfNfa. 
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OiN- 1 ). 
(2.2) 



3 Estimates of the oscillating sum 

The main property of the set is that the Stirling formula holds true on 
this set. On the other hand, the probability distribution P/v(n) is supported 
on this set, up to events with total probability of order 0(N^ 1 ). Hence the 
expectation of the oscillating exponents can be restricted to Bn with the 
same approximation error: 

Z= P N {n)e lSN{n) + 0{N- 1 ). (3.1) 

neB N 

The trapezoidal approximation method is applicable to the sum (3.1) with 
the approximation error (1.3) (see (A.l) for details): 

Z= [ dxf N (x)e l(TN ^ + 0(N~ 1 )+ lN , 

where 

\1n\<1 E {max|Pi 2) (y)|+2max|P^ 

nt-Djv 

By definition of the set D N , we have the following uniform estimates of the 
derivatives 



S^ = ^ + ^ = 0(N^ 2 ), S^ = ^ = 0(N-^), x E D 



'N VN 3 VN 3 
and as is proved in Appendix (see (A. 3) and (A. 4)) 

N 



Pg\x) =P N (x) (log- + 0(x~ 1 )) =P N (x) OiN-W), 
P$\x) = P N (x)((\og^y +x~Hog^+0(x- l \ogx) ) J =P Ar (x)-0(iV- 1 ), 



because N/x = 1 + (iV MogiV) 1 / 2 uniformly on Dn, and hence \og(N/x) = 
0(N~ 1 ' 2 ) and (log(iV/:r)) 2 = O(N^) uniformly in x G D N . 

Since I]„-P/v(n) = 1, it follows that |7jv| = 0(N^ 1 ). Thus we proved that 



_jl 

Z = -j= dxe~ x - , C iS N{ VN X+N) + 0{N -i^ 

\ ^ 2y/W ^ 12N J 

where S N (y/Nx + N) = A(x + y/N) + B(x + y/N) 2 / (2y/N) . Now we expand 
the exponent into the series in xj y/N, and extend the of integration domain 
to the whole real line: 

Z = ^=\ dxe-£(l- WW!v) + ^Wlv)V(2v^) + f N -iy 

y/h^JR V 6^ / 

(3.2) 

The error related to the extension of the Gaussian integral from Dn 
to R can be estimated by an inequality of Komatsu type: 

/ dxx n e~^ < 4(a + l) n ~ (Va 2 + 2 + a) * (3.3) 

J\x\>a \ J 

(see Appendix 4), which implies the required absolute estimate 

1= f dxx n e- ax2/2 < (logiV) n/2 • OiN' 1 ) = OiN" 1 ), 

2TT J\x\eD N a=y/2logN 



27T J\x\eD N 

n = 1,2. 

By using the standard technic, we find elementary Gaussian integral (3.2): 

1 f , _„^2 /OJ _^ ib 



'2-k 



/ dxxe- ax2 / 2+ibx = -^e' b2 / 2a , 
jr. a 6 ' 2 



L [ dxx^ 2 ^ = tb{3a 7 - b2 K ^. 

2tt Jr a 7 ' 2 



Hence, by substituting a = 1 — iB/y/N, b = A + B, we obtain 
Theorem 3.1 

e 2(i-iB/^)(i_j V ^ } )+Q(N- 1 ). (3.4) 
1 — iB/y/N V 6ViV / 
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Figure 3: The picture shows the absolute approximation error E of the sum 
Z by integral (3.2), as a decreasing function of A and B for N = 1000. 

4 Approximation of related series 



In this section we consider the integral approximation of the related series 

n=l 



N ~ N n , ( f An Bn 2 \ 1 . , 



for s = 1,2, which are used to calculate the photo current and its dispersion 
(see @, Chap. 2). Note that Zi can be rewritten as 



N ™N n ( /A(n+1) B(n + 1Y\} 



f+i- 



= Ne^ ^ ■ Z(A + B/N, B, N) = NeVW . z(A, B, N) + 0(1). 
This relation readily implies 

■ A e i(A+B/2)VN {A+B) 2 , (A + B) 3 \ 

Z 1 = Ne^n , e -2(i-iB/VN)(i-i ±flZJ±l ) +0(1). 

'l-iB/VN ^ sVN J 
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Similarly, 

_ ^ N n ( /A(n + 2) B(n + 2) 2 \) 
= Zi + iVV^ (A+B/A ° • Z(4 + 25/iV, 5, TV) + O(iV) = 

■ 2A 

= N 2 e^Z(A,B,N) + 0(N). (4.2) 
For the general case, we can prove the following assertion. 
Theorem 4.1 The sum (4.1) has the following asymptotic expansion 

Z s = N s e^Z(A, B, N) + 0(7V S_1 ), 

where Z(A,B,N) is given by Eq. (3.4). 

Another series which can be estimated with the help of previously devel- 
opped technics is 

W"E^exp{^ + J|)} (4.3) 

Indeed, for the new variable x = (n — N)/y/~N, we have n/y/N = x + \fN, 
and 

y/E = (x + y/N)(l + x/Vn)- 1/2 = y/N + -- + OiN- 1 ). 

Hence, 

A^ + -^L = A ( VN + f - -^=) + B(VN + x)/2, (4.4) 
2^N V 2 8y/NJ 

and by replacing the exponent in Eq. (3.2) by (4.4), we obtain the following 
assertion 

Theorem 4.2 The sum (4.3) has the following asymptopic exapnsion 

iVN(A+B/2) (A+B) 2 . (A + B) 3 \ 

=e ~ sn+MA/AVN) (l-i y ^Zl ) +0(N- 1 ). (4.5) 
V 48^ / 



1 + ZiA/As/N 
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To conclude this section, we consider an estimate of the following multiple 
sum: 

oo jfm+n 



-2N N (.fa 1 m + a 2 n b x m + b x n + 26 3 mn\ 1 , A n . 

z = e L^ exp \i^^ + iTP — )/ (46) 



N m+n r . / aim + a 2 n b\m + b\n + 2b 3 mn\ 

m,n=0 

which was derived in [[| to describe the interaction of the lazer beam in some 
standard state with the quantum oscillator. 

In variables x = (m — N) / VN, y = [n — N) / \^N, this series can be 
approximated by the Gaussian integral 

z = t^H. J w dxdye ^(i - * x+y ^-*) + O(iV-) 

where cx(a, 6) = a\ + a 2 + b 3 + (61 + b 2 )/2 and 

S N (x,y) = -{x 2 + y 2 )/2+ 

+i(a 1 + 61 + b 3 )x + i(a 2 + b 2 + b 3 )y + i^x 2 + b 2 y 2 + 2b 3 xy)j2\fN. 

This integral can easily be evaluated. As a result we obtain the following 
estimate. 

Theorem 4.3 The sum (4.6) has the following asymptopic exapnsion 

e iVN*(a, b )-G N (a,b) , + ^ + ^3 + + ^ + ^3, 

/ = , — 1 - z 7= +C(iv ), 

sjl -i(bx + b 2 )/VN K Q vN J 

where 

, (ai + 6i + 6 3 ) 2 + (a 2 + 6 2 + 6 3 ) 2 , .63(01 + &i + 6 3 )(a 2 + b 2 + 63) 
0Ar(a,6) = === + 1 m . 

2^1 - z(6 a + 6 2 )/ViV ViV 

5 Appendix: Simple a priori estimates 
5.1 Trapezoidal approximation 

Consider the estimate of the trapezoidal approximation of the integral of a smooth 
function f(x). The Taylor expansion implies 

\{f{Xn) + f(Xn+l)) = f(x n+ l) + Kn + K n+ l, 
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where the approximation error does not exceed the value 

Kl < \{x n ~x n+ i) 2 max \f (2) (y)\, 

<Z 2 y€[xn,x i] 



"2 



and for x £ [x n , x n +i], we have 

f(x) = f{x n+h ) + (x - x nH )fW(x nH ) + 7n (x), 

where |7„(x)| < 5(3; — x ra+ i ) 2 max J/e [ a , n jXn+1 ] Hence the error of the trape- 

zoidal approximation can be estimated as follows 

I f(x)dx=-^— -Un + fn+l- K n - K n+ lj - J n (x) dx = 

= Xn+1 ~ Xn (fn + fn + l)+Vn, (Al) 

where /„ = f(x n ) and \i] n \ < \{x n+1 - x n ) 3 max ye[Xn;Xn+l] |/ (2) (y)l- 

5.2 Derivatives of the Poisson distribution 

Consider the estimates of the derivatives of the distribution Pn(x). Note that 
from the Stirling expansion for the Gamma-function 

IX. + 1) - (f ) W(l + ± + ^ + <*»-)) (A.2) 
we have the following estimate for its derivatives: 

r« ( x + 1) = (fog x + + 0(aT 2 )^ T(x + 1) 

and similarly, 

r( 2 )(x + i) = |^io g x + ^ + o(^ 2 )) 2 + Q + o(x- 2 ))|r(x + i) = 

(log 2 x+^ + i + 0(x- 2 logx)lr(x + l). 

L X X J 

Therefore, 

t«(x + i)\ 2 r( 2 )(x + i) n , ( ,„ , 

-0(x _1 logx), (A3) 



r(x + i) / r(x + i) 
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and moreover, 



p (i), v -n d N* f rW(x + i) 

P K N '{x)=e ————— = P N (x)[ log N 



dxT(X + l) JVV J \ & r(x+l) 
P w (x)(log^-i- + 0(x- 2 )Y (A4) 



Similarly, for the second derivative we have 

d 2 N x 



P^(x) = e 



dx 2 T(X + 1) 

2 r( 2 )(x + i) /rW^ + m 2 * 



r / A" i ^. 2 a z r^ i + i /rw i + i vi 

= Fw(x){ flog ^ - i + Ott- 2 )J + 0(i-' log i) J . 
Finally we obtain, 

£ i#<»)- E P«(»){(lo g ~) 2 -M^ + 0(n -Mo g „)}. 

~ _ / 2 log AT / 2 log iV ~ i 9 /V I 

where Bn = [Ne V jv 7 Nev n ]. Note that for n G -Bat, we have | log — | < 
lo f} N , where log N = 0(1). By similar reasons, log ^" — - + 0(n~ 1 logn) = O^N^ 1 ) 
uniformly on this set; in the sense of asymptotical expansions, Bjy is equivalent to 
the previously defined set Bjq. Hence, 

J2 Pg\n) = OiN- 1 ) J2 Pn(u) = OiN- 1 ). (A5) 

5.3 Corrections to the Gaussian distribution 

From the Stirling expansion (A. 2) we have 

p N ( n ) = e -N+n-n\ og % 1 

y ' V2^N v ^7]V(l + (12n)- 1 + 0(n- 2 ))• 

Set A" = x\fN + N; then for x G Dn we obtain 



-JV + n - n log — = -N + xVN + N - (xVn + N) log(l + x/VN) 
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nr? i /T7 f x x 2 x 3 x A ^( x 5 \ 
xVN - (xVN + N)\ — = - — + — = + — ? + O -= 



x 2 X 3 x 4 / X 5 



- 2 + ^N-Tm + {7m)- (A6) 

This expansion contains corrections to the Gaussian exponent as in (1.1). Similarly, 



2VN 




V^(l + (I*.)- + O(n-)) = ^1 + -^(l + J^(i - + 0(w -')) 

-}j + 0(x/iV 3/2 ) +0{N- 1 ij = 



+ 0{N- 1 ) + 0(x/N^ 2 ). (A.7) 



Equations (A.6)-(A.7) imply expansion (1.1). 

5.4 Komatsu inequalities 

The inequalities of the Komatsu type can easily be justified by the following con- 
struction. For a positive absolutely integrable function /, we consider the integral 
Sx° f(y)dy = F(x), F(oo) = 0. If there exists some smooth positive decreasing 
function g{x) such that g{oo) = and 

then F(x) < g(x). Indeed, the inequality (A. 8) for derivatives implies that g 
decreases at infinity more rapidly then F. Since they coincide at infinity, this 
proves that F(x) < g{x). In the Komatsu inequality and in our generalization 
(3.3), the right-hand sides satisfy (A. 8). For the left-hand sides, the proof is 
similar; it uses the characteristic inequality opposite to (A. 8). 



References 

[1] W. Feller, An introduction to probability theory and its applications, 
1, John Wiley & Sons, 1977. 

[2] S. Bose, K. Jacobs, P. L. Knight, "Preparation of nonclassical states in 
a cavity with a moving mirror", Phys. Rev., A 65, N5, 4175-4186. 



14 



[3] C. Brif, A. Mann, "Quantum statistical properties of the radiation field 
in a cavity with a moving mirror" , Preprint LANL, |quant-ph 9905040| , 
9p. 

[4] A. M. Chebotarev, Lectures on Quantum Probability, Aportaciones 
Matematicas, Sociedad Matematica Mexicana, Textos Avansado, N14, 
Mexico, 2000, 294 p. 

[5] K. Ito and H. P. McKean, Diffusion processes and their sample paths, 
Springer, 1965. 



15 



